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SUMMARY 

A major life-limiting phenomenon of electric thrusters Is the sputter 
erosion of discharge chamber components. Thrusters for space propulsion are 
required to operate for extended periods of time, usually In excess of. 

10 000 hr. Lengthy and very costly life-tests In high-vacuum facilities have 
been required In the past to determine the erosion rates of thruster compo- 
nents. It Is, therefore, highly desirable to develop alternative methods for 
determining erosion rates which can be performed In relatively short periods 
of time at considerably lower costs. 

This paper presents an attempt to relate optical emission Intensity from 
an Ion bombarded surface (screen grid) to the sputtering rate of that sur- 
face. The model used a kinetic steady-state (KSS) approach, balancing the 
rates of population and depopulation of ten low-lying excited states of the 
sputtered molybdenum atom (Mol) with those of the ground state to relate the 
spectral Intensities of the various transitions of the Mol to the population 
densities. Once this Is accomplished, the population density can be related 
to the sputtering rate of the target. Radiative and colllslonal modes of 
excitation and decay were considered. Since actual data has not been pub- 
lished for Mol excitation and decay rate constants, semlemplrlcal equations 
were used. 

The calculated sputtering rate and Intensity was compared to the measured v 
emission Intensity of the excited sputtered Mol and sputtering rates of the 8- 
and 30-cm Ion thrusters using semlemplrlcal equations for the electron density, 
temperature of the discharge plasma and sputtering yields of a Hg -> Mo system. 


INTRODUCTION 

A major life-limiting phenomenon of electric thrusters Is the sputter 
erosion of discharge chamber components. Thrusters for space propulsion are 
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required to operate for extended periods of time, usually In excess of 
10 000 hr. Lengthy and very costly life tests In high vacuum facilities have 
been required In the past to determine the erosion rates of thruster compo- 
nents. Alternative methods for determining erosion rates which can be 
performed In relatively short periods of time and at considerably lower costs 
are, therefore, highly desirable. 

This paper presents an attempt to relate the optical emission Intensity 
from an Ion bombarded surface (screen grid) to the sputtering rate of that 
surface. Several methods have been tried and tested. One consists of attach- 
ing a "badge" consisting of alternate layers of very thin metal layers to the 
sputtered surface (ref. 1). Such a technique has not always yielded results 
which agreed with long term life tests (ref. 2). Another technique would be 
to apply absorption spectroscopy to the life-time determination. Emission 
spectroscopy was applied by Wehner about 20 yr ago In an attempt to measure 
relative sputtering rates of various target materials by Inert gas Ions at low 
energies (ref. 3). The Ion thruster plasma, however, represents a much more 
complicated system than the experimental arrangement of Wehner. 

First, the thruster discharge plasma was characterized In regard to the 
equilibrium conditions and optical depth of the plasma. The lack of thermal 
equilibrium complicates the understanding of the plasma as well as the formu- 
lation of the model. The determination of the optical depth of the plasma 
determines If any radiation Is absorbed within the plasma. 

This work uses a steady state kinetic model to Infer the population dynam 
Ics of the sputtered molybdenum "Impurity" atoms In an Ion thruster plasma. 
This approach requires a knowledge of basic plasma parameters as well as spec- 
troscopic emission data of the Important molybdenum atomic transitions. A 
total molybdenum atom population can be found, which can be related to the 
sputtering rate of the molybdenum surface, using known relationships. The 
model formulated was compared to measured relative spectroscopic Intensities 
from an 8- and 30-cm Hg Ion thrusters. For a more complete treatment of this 
study see reference 3. 


CHARACTERIZATION OF ION THRUSTER PLASMAS 

The parameters, which are of primary Importance In characterizing the 
properties of a plasma are Its electron temperature and electron density. In 
many plasmas, the energy distribution of the electrons Is sufficiently close 
to a Maxwellian form so that a meaningful temperature can be assigned to the 
electrons; this temperature may not be equivalent to the temperature of either 
the atoms or the Ions In the plasma, particularly In some laboratory plasmas 
(ref. 4). The range of electron densities In plasmas varies from below 10 4 
to above 10 17 electrons/cc . Electron temperatures, typically expressed as 
energies, range from 0.01 to over 100 000 eV. Such a wide variety In these 
basic parameters Implies a wide diversity of plasma characteristics. For the 
purpose of characterizing plasma behavior the colllslonal-radlatlve model 
covers virtually all known plasmas. There exist special cases for which sim- 
plifications In the required calculations can be made. Two of these special 
cases are the corona and the local thermodynamic equilibrium models. 

Additionally, the optical thickness of a plasma must be considered. An 
optically thin plasma Is easier to treat for both the colllslonal-radlatlve 
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and the corona models. Both of these models are based on the assumption that 
the plasma Is optically thin, particularly with respect to Its resonance lines. 
Where this assumption Is not valid, corrections must be added to account for 
the plasma's absorption (ref. 5). 


The Local Theromodynamlc Equilibrium Model 

The simplest of the models from a calculatlonal point of view Is the local 
thermodynamic equilibrium (LTE) model. The LTE model assumes that colllslonal 
processes are dominant In the plasma and that the particle collision rates are 
high enough so that the system can Instantaneously adjust to any change In 
plasma conditions. A system In local thermodynamic equilibrium Is governed by 
the same mathematical expressions as systems which are In total thermodynamic 
equilibrium although, an LTE system does not meet all of the criteria required 
of a system In total thermodynamic equilibrium. In an LTE system, population 
densities need not be equal at different points In the plasma. The plasma 
also need not have a single overall temperature. At any point In the plasma, 
a local temperature Is meaningful and the density and chemical composition of 
the plasma at this point. In conjunction with this temperature, can be used to 
define the population of the various energy levels In that region by means of 
a thermodynamic equilibrium approximation. 

If colllslonal processes are to dominate radiative ones In a plasma, suf- 
ficient number density must be present so that frequent colllslonal excitation 
and decay can occur. If an optically thin plasma Is assumed, then radiative 
excitation can be Ignored. The criteria for the existence of local thermo- 
dynamic equilibrium would then be that the colllslonal decay rate must be sig- 
nificantly greater than the radiative decay rate. 

A wide variety of possible colliding species exist In an Ion thruster 
plasma. It should be noted that most of the species are present to a much 
smaller extent than the two major components of the plasma, propellant atoms 
and electrons. These two species are present In concentrations which are 
usually within an order of magnitude of each other and which exceed by at 
least two orders of magnitude the populations of any of the other species. 
Electrons are much more efficient In producing excitation and decay than are 
atoms, for reasons to be discussed In a later section. Therefore, It Is rea- 
sonable to assume that colllslonal effects on excitation and decay are solely 
governed by collisions between plasma electrons and plasma atoms. It Is for 
this reason that electron density and energy are Important plasma parameters. 

Sputtered molybdenum atoms In low-lying excited states, are also of 
Importance In the determination of both the rate of colllslonal decay and In 
the rate of radiative decay; both rates are proportional to the population 
densities of the excited molybdenum atoms. The colllslonal decay rate Is the 
product of the electron density of the plasma, the density of molybdenum atoms 
In the excited state being considered and a colllslonal deexcitation rate con- 
stant, yielding a transition rate per unit time and unit volume. The corres- 
ponding radiative decay rate Is the product of the density of the excited state 
of the molybdenum atom and Its Einstein transition probability. A criterion 
often used to determine the validity of the LTE assumption Is that, for LTE to 
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exist, the colllsional decay rate must exceed the radiative decay rate by at 
least a factor of 10: 


n e n j R CD - 10 ' n ^ ^ A 


j X 

1<J 


or, equivalent to this, ' 


"co 2- 10 k rd 


(l) 


where n e Is the electron density In electrons/cc, nj Is the density of the 
j (upper) state of the molybdenum atom, In atoms/cc, RcD Is the colllsional 
deexcitation rate constant, In transitions s -1 electron -1 atom -1 , and Aj^ 

Is the Einstein transition probability, In transitions atom -1 s -1 (often 
expressed as s -1 ). The sum of these transition probabilities below state j 
represents the radiative decay of state j to all lower states. When non- 
resonance transitions are Ignored, this summation reduces to a single transi- 
tion probability Aj^, where the state 1 Is taken to be the ground state of 
the atom. «cd represents the rate of colllsional decay and krq represents 
the rate of radiative decay, both In units of transltlon/s. When the above 
Inequality Is satisfied, a local thermodynamic equilibrium approach Is 
justified. 


It can be seen from equation (1), that plasmas with very high electron 
densities are more likely to be In local thermodynamic equilibrium than low 
electron density plasmas (ref. 5). 


The Corona Model 

Another Important model of plasma behavior Is the corona model. This 
model was proposed to explain some of the features of the solar corona; It Is 
equally applicable to low-density laboratory plasmas. 

The corona model Is based on the fact that, In a low-density plasma, 
radiative decay will predominate colllsional decay. This occurs because a 
low-density plasma cannot support a rate of collision high enough to compete 
with radiative decay. The assumption of optical transparency can usually be 
made for the corona model. The criterion for coronal equilibrium Is therefore 
that the colllsional excitation rate must balance the radiative decay rate. 
Equation (2) states this mathematically: 

n e n 1 R CE = n j ^ A j1 (2) 

1<j 


or, equivalently, 


k CE = k RD 
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where Is the density of the molybdenum atom lower state In atom/s, Rr^ 

Is the collislonal excitation rate constant. In transitions s _1 electron-' 
atom -1 , and *ce represents the rate of collislonal excitation. Other terms' 
are as defined In equation (1). 


The assumptions of the coronal model break down as the density of a plasma 
goes up. When collislonal decay becomes significant the balance between col- 
llslonal excitation and radiative decay Is lost. A criterion similar to that 
given to evaluate local thermodynamic equilibrium can be developed for the 
coronal model. When the collislonal decay rate exceeds 10 percent of the 
radiative decay rate, the corona model becomes Inappropriate. 

n j ^ A j1 - 10 {n e n j R CD ) (3) 

1<J 

or, equivalently, 


K 


R0 - 


K 


CD 


When the above Inequality Is satisfied, a corona equilibrium approximation Is 
justified. Corona equilibrium Is favored by low electron densities (ref. 24). 
As will be shown later, over the range of electron temperatures and densities 
of an Ion thruster considered, radiative decay always exceeds collislonal 
decay. This Implies that the corona model may also be appropriate for an Ion 
thruster. 


The Colllslonal-Radlatlve Model 

For those situations where neither the local thermodynamic equilibrium 
model or the corona model are appropriate, the colllslonal-radlatlve model can 
be used. This model accounts for the effects of both collislonal and radi- 
ative decay on the upper levels of the system. To some degree. It bridges the 
gap between the LTE model, where collislonal decay dominates the decay process 
and the corona model, where radiative decay Is predominant from all excited 
states. The colllslonal-radlatlve model Is appropriate for those plasmas which 
have an Intermediate electron density, such that neither radiative decay, nor 
collislonal decay can be neglected (ref. 5). This method will be applied 
herein. 


Optical Thickness, Absorptivity 

It was stated earlier that the opacity of a plasma to Its various spec- 
tral lines could complicate the Implementation of a plasma model. In particu- 
lar, radiation trapping In the corona and colllslonal-radlatlve models can 
Increase the calculatlonal overhead considerably. Some method Is needed to 
quantitatively evaluate the optical thickness of an Ion thruster plasma at any 
given frequency or wavelength of radiation. A commonly used approximation Is 
the method described by Grlem (ref. 6). This method makes several assumptions 
about broadening mechanisms In plasmas as well as values for "typical" oscil- 
lator strengths In plasmas, which are not appropriate to the Ion thruster sys- 
tem. An alternate method of evaluating optical thickness, detailed by Alkemade 
(ref. 7) and without the limiting assumptions of Grlem will be discussed. 

5 



For a uniform plasma, a quantity called the absorptivity can be defined as 

L(v q ) = 2.65X10" 2 ^ f 1J /a« eff (4) 

where L(v 0 ) Is the absorptivity of the line at line center, In cm -1 , n^ Is 
the number density of the lower state of the atom In atoms/cc, fM is the 
absorption oscillator strength for the transition under consideration and 
Av e ff Is the Doppler width at half-intensity In centimeter (ref. 7). Besides 
allowing pertinent oscillator strengths to be used, this approach permits a 
choice of methods to estimate the spectral line broadening parameter, Av e ff. 

In the case of the Ion thruster, Doppler broadening will be present, but In a 
nonstandard form. Normal Doppler broadening can be estimated using 
equation (5a), 

5w q = 2(u q /c)V( 2 In 2) kT/m (5a) 

where 6 vq Is the normal Doppler spectral halfwidth In s _1 , k Is the 
Boltzmann constant, T Is the atom temperature In K, m Is the mass of the 
species of Interest, In g, v 0 Is the frequency at the line center In s' 1 
and c Is the speed of light In cm/sec (ref. 8). 

Normally, the result from equation (5a) would be used In equation (4) to 
calculate absorptivity. However, Doppler broadening Is expected to be atypical 
for molybdenum atoms In an Ion thruster (ref. 3). The following equation 
approximates the effect of this difference. 

Aw eff*2 5v D (5b) 

The result from equation (5b) Is used In equation (4), yielding the values 
shown In table I. 

To evaluate optical thickness using absorptivity, the absorptivity Is 
multiplied by 8., the path length of the plasma along the optical axis. If 
this quantity Is significantly below one, the plasma can be considered opti- 
cally thin. The recommended upper limit for L(v 0 ) 8. Is between 0.1 and 0.2. 
This allows for the error present In the values used for these calculations 
(ref. 9). Typical values for 8 In the Ion thruster are on the order of 

10.0 cm. The data In table I shows that, for molybdenum kinetic energies above 

1.0 eV, the plasma Is transparent to the radiation of the ten transitions dealt 
with by this model. Evidence will be presented later showing that molybdenum 
energies will exceed this 1.0 eV limit and therefore, the plasma can be con- 
sidered optically thin. 


FORMULATION OF KINETIC STEADY-STATE (KSS) MODEL 

The goal of this work Is to relate measured Intensities of molybdenum 
spectral lines, to the rate of sputter damage to the molybdenum screen grid of 
the Ion thruster. It Is desirable to keep the model as simple as possible, to 
avoid the need for large amounts of computer calculations. It Is also Impor- 
tant to keep the model as general as possible, to allow It to be adapted easily 
for use on other types of electron bombardment Ion thrusters, with different 
propellants, such as xenon, or with different magnetic field configurations. 
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The model assumes that the molybdenum Ionization rate Is so small that It 
can be neglected. A number of molybdenum atom transitions have been seen In 
the spectrum obtained from a typical 30-cm thruster. The Inclusion of all 
these transitions Is not necessary If the ground state of the molybdenum atom 
predominates. If more than 90 percent of the molybdenum atoms are In the 
ground state, the use of a subset of the observed molybdenum transitions 
should be quite satisfactory. 

The model considers the molybdenum ground state and ten low-lying excited 
states, which are known to have resonance transitions In flames (ref. 11). 

Only resonance transitions are considered herein. This Is done for a number 
of reasons. The resonance transitions should be significantly more Intense 
than the nonresonance transitions In a relatively low temperature plasma 
because of the much lower population densities of the excited states, which 
are the lower states In all nonresonance transitions. Spectroscopic data, 
such as oscillator strengths and Einstein transition probabilities, for molyb- 
denum atom nonresonance transitions tend to be either of poor quality or 
unknown (ref. 12). In addition, most of the nonresonance transitions In the 
eleven state system under consideration will be In the IR region, where accu- 
rate Intensity measurements are very difficult. For these reasons, nonreso- 
nance transitions were omitted from the model. 

Two major steps are required to achieve the goal of relating the observed 
Intensity of molybdenum atom transitions to the screen grid sputtering rate. 
The first step Is to relate the spectral Intensities of the various transi- 
tions of the molybdenum atom to the population densities of the various states 
of the atom. The second step Is to relate this population density data to the 
screen grid sputter rate. 


Steady-State Theory 

The task of calculating the population densities of the various energy 
levels of a species Is relatively easy for systems which are In thermodynamic 
equilibrium or which satisfy the criteria for local thermodynamic equilib- 
rium. Such a solution relies on the fact that the level populations of a sys- 
tem In or near complete thermodynamic equilibrium behave In a statistically 
predictable manner, depending only on the temperature of the system , the 
statistical weights of the levels and the energy differences between these 
levels. All the particles In a system which are In complete thermodynamic 
equilibrium share a common temperature, that Is, each of the species In the 
system has a Boltzmann distribution of particle energies which corresponds to 
one common temperature (ref. 13). Systems which are In local thermodynamic 
equilibrium share a common temperature In a given region of the plasma; a 
different region of the plasma may have a different characteristic tempera- 
ture, but all the particles In that region share that same temperature 
(ref. 5). Systems which are not near thermodynamic equilibrium are more 
complicated to treat since they do not behave In the normal statistically 
predictable manner. For these systems, Individual rates for each of the 
processes occurring In the plasma must be considered (ref. 14). Such rates 
have not, In general, been measured for most species; usually empirical or 
semlemplrlcal data Is used to estimate the rates of these processes. 
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A system such as the Ion thruster plasma, governed by a number of pro- 
cesses, should eventually reach a point where the rates of all processes which 
populate any given level will be balanced by the rates of all processes which 
depopulate that level. When such a situation occurs, the population of each 
level Is In a steady state condition. As long as no parameter Important to 
any of the processes changes, the population of each level will remain con- 
stant. The steady-state approximation allows a system of rate equations to be 
established. When this system Is solved, values for the level populations are 
generated and In theory, the first major step toward the stated goal of this 
work Is achieved (ref. 15). 


Excitation-Decay Rate Constants 

A number of simplifying assumptions were made previously about the nature 
of an Ion thruster plasma. One of these assumptions Is that, since molybdenum 
Ion lines have not been observed In Ion thruster spectra. Ionization can be 
considered negligible. In such a case, only atomic excitation and deexcita- 
tion processes need to be considered. This eliminates from consideration the 
large numbers of Interactions In which Ions can be Involved In a plasma dls 
charge (ref. 16). This leaves two predominant modes of excitation and decay. 

Excitation 

1. Radiative excitation 
2. Colllslonal excitation 

Decay 

1. Radiative decay 
2. Colllslonal decay 

Radiative excitation occurs when a photon Is absorbed by an atom of one 
electronic energy state, yielding an atom at a higher energy state. This 
process can occur In plasmas and flames with a very high spectral radiation 
density, such as those which are optically dense or which are Illuminated by 
an external light source. Photons which are absorbed by an atom must be very 
close In energy to the energy required for one of the electron transitions In 
the atom. A plasma In which photon absorption Is occurring will often be 
optically transparent over most of Its spectrum, but In the vicinity of some 
of Its resonance lines. It can be optically thick (ref. 17). Any assumptions 
of optical transparency should be checked using the methods detailed previ- 
ously. Results from table I show that the Ion thruster plasma Is optically 
transparent to all the molybdenum lines used In this study. It should be 
noted that the plasma will appear optically thick to the mercury resonance 
line listed In table I. Since mercury atom Intensities do not enter Into the 
model, this result does not cause a problem. This value for mercury was 
Included In the table only for comparison. Since the molybdenum atom reso- 
nance lines under consideration satisfy the criteria for determining optical 
transparency by a sufficient margin, radiative excitation should be negligible. 

In light of the previous argument against radiative excitation, colll- 
slonal excitation would appear to be the only excitation mechanism of Impor- 
tance In the main discharge plasma. The various species with which a 
molybdenum atom may collide were briefly considered previously. Such a colli- 
sion might Involve atoms, Ions or electrons. The atomic species present In 
the thruster are mercury, molybdenum and contaminant atoms. Mercury Ions will 
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also be present, at a density roughly one to two orders of magnitude less than 
mercury atoms (ref. 18). Most of these Ions will be singly Ionized, with a 
small, but significant number of them being doubly Ionized. There may be also 
trace amounts of molybdenum and contaminant Ions In the plasma; such species 
should be present In negligible amounts, compared to mercury Ions. Electrons 
are also present In the discharge; they should be present In the plasma In 
numbers roughly equivalent to the population density of mercury Ions In the 
plasma (ref. 21). A list of these species, In descending order of their 
expected population density Is given here: 

Hg atoms, electrons > Hg Ions > Mo atoms > 
contaminant atoms > Mo and contaminant Ions. 

It Is unlikely that any species less common that mercury Ions will cause 
significant colllslonal excitation. The molybdenum atom number density Is 
expected to be between three and five orders of magnitude less than mercury 
atoms (ref. 19). In general, atom-atom and atom-ion collisions are not very 
effective In exciting electronic transitions. Atomic excitation by collision 
Is best thought of In terms of a collision with one of Its bound electrons. 
Energy transfer between particles of such widely differing masses Is very 
Inefficient. Transfer of the proper quantity of energy Is very unlikely (ref. 
20 ). 


As any molybdenum atom which does not collide with a mercury atom (Ignor- 
ing the other far less populous atoms and Ions) will collide with a thruster 
surface and stick. It appears that not only will atom-atom (or atom-ion) col- 
lisions not be very productive In producing excitation, they will also not be 
particularly common. 

Similar reasoning can be applied to electron-atom collisions to show that 
these collisions should be efficient In producing excitation. A collision of 
a free electron with a bound electron can result In the transfer of a large 
fraction of the energy of the free electron to the bound electron (ref. 21). 
The Maxwellian population of electrons have an energy on the order of several 
electron volts. Since the energy differences between the levels being con- 
sidered In this model are also several electron volts, there will be a rela- 
tively large number of electrons of the proper energy available to excite such 
transitions. Electron-atom collisions are therefore, considered to be the 
only Important source of excitation of molybdenum atoms In the plasma. The 
radiative decay rate Is dependent only on the Einstein transition probability, 
Aj^, and the population of the upper state of the atom (ref. 15). 

Colllslonal decay Is the Inverse process of colllslonal excitation. The 
same reasons which were used to justify the exclusive Importance of electron- 
atom colllslonal excitation apply to colllslonal decay. Efficient energy 
transfer and relatively high electron densities make electrons the only col- 
lision partners capable of causing significant decay of excited states. Col- 
llslonal decay might be expected to play a lesser role than radlatlonal decay 
In many situations; the low number density of the excited states gives the 
electrons less opportunity to exchange energy In collisions. This Is espe- 
cially true If the electron density of the plasma Is relatively low. 

Since actual data has not been published for molybdenum atom excitation 
rate constants, semlempl rlcal equations were used to estimate the values of 
these constants. These equations are 
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(1) colllslonal excitation 


-5 ^11 / x \ 

R - 1 .58x10 13 ** exp (- *-V (6) 

xT e ' e/ 

where Rqe Is the colllslonal excitation rate constant, f^j Is the oscilla- 
tor strength for absorption, a unitless term, <G| C a > > 1s the thermally averaged 
Gaunt factor, a unitless term x Is the transition energy In electronvolts , T e 
Is the electron temperature, also In electronvolts and n e Is the electron 
density, In electrons/cm 3 (ref. 15). 

(2) colllslonal decay 


R cd = 1 .58x10 


-s/M lu/V 


XT 


1/2 


(7) 


where Rco Is the rate constant for colllslonal decay, g^ Is the statistical 
weight of the lower electronic state and gj Is the statistical weight of the 
upper electronic state (ref. 15). 

(3) radiative decay 

r r D = ( 8 ) 

where Rrq Is the rate constant for radiative decay and Aj^ Is the Einstein 
transition probability (ref. 15). 

These rate constants permit a steady-state model to be formulated, allow- 
ing an estimate of the'populatlons of the eleven states Included In this model. 


The Kinetic Steady-State (KSS) Model 

Although eleven states of the molybdenum atom are used In estimating the 
molybdenum atom population In the plasma, a four level, atomic system will be 
used to simplify the presentation of equations In this section. Figure 1 
shows this arbitrary system, composed of a ground state and three low-lying 
excited states. The states have been numbered so that state 1 corresponds to 
the ground state and state 4 corresponds to the highest excited state con- 
sidered. All of the following equations can be easily extended to eleven or 
more levels. 

The steady-state assumption requires that after a sufficient time, the 
rates of change of the populations of all levels of a system will balance, so 
that no net changes will occur, In the absence of a perturbation which affects 
the rate of one of the Important processes. This can be stated mathematically 
for the simplified system as 
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( 9 ) 


ll n l 

+ 

R 1 2 n 2 

+ 

R 1 3 n 3 

+ 

R 14 n 4 = 0 

21 n l 

- 

$ 22^2 

+ 

R 23 n 3 

+ 

R 24 n, = 0 

31 n l 

+ 

R 32 n 2 

- 

R 33 n 3 


R 34 n 4 ' 0 

41 n l 

+ 

R 42 n 2 

+ 

R 43 n 3 

- 

R 44N * 0 


where R^j Is the overall rate constant for populating state 1 from state 
j. When 1 equals j, the rate constant Is a depopulation rate constant; 
when 1 Is not equal to j, the rate constant Is a population rate constant. 
Note that the subscript order convention for the rate constants and rates Is 
the opposite of the convention used for spectroscopic terms, such as f^j 
and Aji- This allows the standard mathematical convention of row, column 
subscripts for matrices to be preserved. Only the rates and rate constants 
share this convention. The nj In this equation Is the population number 
density for state j. These are the quantities to be solved for by this 
model. Each equation balances the rate of population of a level 1 (ER^nj) 
with the rate of depopulation of that level (-Rjjnj); the sum of these two 
terms must, under steady-state conditions, equal zero (ref. 15). 


A specific example of the calculation of the overall rate constants for 
both population and depopulation can be Instructive. The subscripts used for 
the partial rate constants are In the form "process, transition". As an exam- 
ple, R CE 21 the partial rate constant for the process for colllslonal 
excitation for the transition which populates level 2 from level 1. Here are 
three examples of the calculations of populatlon/depopulatlon rate constants. 


(a) R ?1 = R ce 2 -| (upward transition, population 

(b) R. „ = R„ n + R ori , 0 (downward transition, 

12 iu,i2 hu, 12 level 2 ) 

(c) R 22 = R 12 + R 32 + R 42 (° vera ^ depopulation 


of level 2 from level 1 ) 
population of level 1 from 

rate constant) 


= (R CD,12 + R RD,12 ) * R CE,32 + R CE,42 


( 10 ) 


The overall depopulation rate constant for level 2 (eq. (10c)) warrants 
an explanation. The depopulation rate constant for level 2 Is equal to the 
sum of all the overall rate constants, by which level 2 populates all other 
levels; the first equation represents this equality. The second equation 
expresses the overall depopulation rate constant In terms of partial rate con- 
stants. Since R-|2 Is a downward transition; the first two partial rate 
constants express the effects of both colllslonal and radiative decay on that 
transition. R 32 and R 42 are upward transitions, the overall rate constant 
for each of these transitions Is equal to the partial rate constant for colll- 
slonal excitation for that process. Each overall rate constant Is created In 
this way to form the rate constants In equation set (9). 


11 



Equation set (9) can be rewritten In matrix form. 


h ii 

R 12 

R 13 

R h\ 


fa\ 

1 

f °\ 

R 21 

- R 22 

R 2 3 

R 2 4 


", 


0 

R 31 

R 3 2 

- R 33 

R 34 , 


n 3 


0 

\ R 41 

R 42 

R 43 

-*u) 


w 


v° y 


This matrix equation Is composed of a rate constant matrix, the population 
vector and a zero vector. The use of standard methods to solve this equation 
typically yields a result of zero for all the members of the population vector. 
This trivial solution occurs due to the zero vector. Perturbing one of the 
elements of the zero vector slightly yields nonzero results for the members of 
the population vector. Unfortunately, these values are dependent on the value 
used to perturb the zero vector. If the value used for the perturbation Is 
small compared to the actual populatlon/depopulatlon rates occurring In the 
plasma, the results obtained for the nj values will be correct In a 
relative sense. It Is conventional to normalize these relative population 
values (the nj vector) so that they sum to one. 


Atomic theory predicts that an absolute population density of an upper 
electronic level can be predicted using the absolute Intensity of a transition 
Involving that upper state 


4irl0 7 1^ 
n j = h „ 0 lAj, 


( 12 ) 


where nj Is the absolute population of the upper state j In atom/cm 3 , Ij^ 

Is the absolute Intensity for the transition j ■* 1 In W cm -2 ster-1, hv 0 
is the energy of the transition In ergs, 8. Is the path length over which 
emission occurs, In cm, and Aj^ Is the Einstein transition probability In 
transition atom-1 S -1 (ref. 39). A single transition between any two of the 
eleven states considered In this model allows the relative population values to 
be converted Into absolute population values using equation (12). A nj value 
Is generated for the observed transition j -» 1 . The ratio between an absolute 
population density, nj, and a desired absolute population density, n^, Is 
equal to the corresponding ratio between their relative population densities, 



(13) 


therefore. 


n 


k = 


Hi 


(14) 


If the desired n^ value Is the total molybdenum atom population density, 
nj, equation (14) becomes 
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since n|< = ny = 1 « If the relative population densities have been normalized 
to one. 

It would be advantageous to use several absolute Intensity values, allow- 
ing several predictions of the conversion factor needed to change nk to nk- 
General agreement of the resulting conversion factors Indicates the model Is 
correctly predicting the behavior of the physical system. 

Tables II (a) to (j) present the calculated nj values from program 
KSS. Each table contains the data for one value of electron temperature and 
several values of electron density. The data are presented In an exponential 
format for each of the terms, 7.15e-7 stands for 7 . 1 5x1 0 -7 . Since these are 
relative population values, they are dimensionless. As can be seen from these 
tables, the percentage of atoms In the ground state never falls below 99.8 per- 
cent for any combination of electron temperature and density; this value lies 
well above the limit of 90 percent mentioned previously. The use of eleven 
levels should thus adequately model the spectral behavior of molybdenum In the 
Ion thruster plasma. 

For electron temperatures and densities which are within the bounds of 
these tables, but are not explicitly tabulated, linear Interpolation can be 
used to find an approximate value for nj . 


APPLICATION OF KSS MODEL TO EVALUATION OF SCREEN GRID LIFETIME 

As has been stated, the goal of this study Is to find a method of corre- 
lating real-time spectroscopic data from Ion thruster emission to the sputter- 
ing (lifetime) of Its screen grid. The formulated KSS model has provided a 
means of calculating nj , the fraction of excited atoms of sputtered molybde- 
num as a function of electron temperature and density. In this section, ni 
will be related to the sputtering rate of the target as a function of the Ion 
thruster operating parameters. Before the theory can be compared to measured 
emission Intensity data from an Ion thruster, the sputtering rate of the screen 
grid has to be determined by calculated or measured sputtering yields as well 
as the plasma properties of the 8- and 30-cm thruster discharge chamber defined 
as a function of the thruster operating parameters. In addition, this discus- 
sion should also lead to a better understanding of the problems associated In 
evaluating thruster lifetime In general. 


Sputtering Rate Calculations 

The sputtering of the upstream side of the Ion thruster has been deter 
mined to be one of the main life-limiting phenomena of Ion thrusters (refs. ?2 
to 24). The sputtering rate of an Ion thruster screen grid bombarded by Ions 
can be calculated from equation (18) (refs. 24 and 25): 

R = [2S f ( E) + (J*VJ*)S^(E)] 

2qnF<p s A B ^ + 


( 18 ) 



where: 


J D beam current 

0 

q electronic charge 

n density of target material 

F beam flatness parameter 

<p s screen grid effective transparency 

A d beam area 

D 

S f (E),S ++ (E) sputter yields of singly and doubly charged Ions 

ratio of doubly to singly charged Ion current density 
E Ion energy 

Beatle (ref. 25) has derived the following empirical expressions for E and 
<p s for a 30-cm mercury Ion thruster. 

E = 1 .28 V D - 9.38 (volts) (19) 

(P S = 0.78 + 0.084 V T - 0.062 J B (20) 

where: 

Vq discharge voltage, V 

V] total accelerating voltage, kV 

These expressions are tfot very sensitive to the operating parameters and, 
therefore, the assumptions will be made for the 8-cm thruster that the Ion 
energy Is equal to the discharge voltage and the grid transparency Is equal to 
the physical open area of the screen grid. 

The sputtering rate R can be also related to the sputtered atom total 
density, nj, by: 

R = n T V Ho 6 (21) 

where 


v Mo average sputtered atom velocity 
6 cross sectional area of the discharge 

By combining equations (12), (15), and (23) the sputtering rate can be finally 
related to the emission Intensity of a particular transition as: 


R 


4 « x1 ° I il V Ho 6 

Vj1 h V 


( 22 ) 
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The relative sputtering rates at two thruster operating points thus becomes: 


»1 _ Nib); 

R 2 l) 2 ( nj )-| 


(23) 


Vfjj 0 In equation (21) represents the sputtered molybdenum atom velocity 
term. This term Is expressed as an energy because that Is a customary method 
of reporting this value In the literature. Finding suitable absolute values 
for this term Is not a trivial matter. Various researchers (refs. 26 to 32) 
have determined that this parameter ranges from a few to several tens of elec- 
tronvolts for most target-ion combinations at the relatively low Impact ener- 
gies characteristic of Ion thrusters. A value for the most probable ejection 
energy of the sputtered target material can be derived from collision cascade 
theory (ref. 32). This value, Independent of Incident Ion energy, Is roughly 
one-half the surface binding energy of the target material (ref. 28). Since 
the surface binding energy of molybdenum Is 6.82 eV (ref. 33), the most proba- 
ble energy for sputtered molybdenum atoms would be 3.4 eV. The energy distri- 
bution of sputtered metal atoms tends to obey the following relationship 


P (E) = 1 (24) 

(e + u r 

where p(E) Is the differential flux of sputtered particles with an energy of 
t, and U Is the surface binding energy of the metal, In consistent units of 
energy. This distribution has a pronounced high energy tall; the average 
energy of sputtered atoms should therefore, be significantly higher than the 
most probable energy. 

The sputtered atom velocity, Vfu| 0 , has been assumed to be Independent of 
Incident Ion energy. Should direct measurements of this parameter be reported 
In literature, they can be easily Incorporated In the model. 


To evaluate the spectroscopic model formulated In this study only rela- 
tive sputtering rates are necessary. Therefore, for the same thruster config- 
uration, using equation (18), the relative sputtering rates become: 


where 


Rl Vs 2 ( S < E »1 

R 2 = ^sl^^ 


(25) 



Considerable amount of sputtering data has been accumulated over the years 
(refs. 34 and 35). However, reliable data, especially at Ion energies near 
the sputtering threshold (region of operation of an Ion thruster) Is still 
rather scarce. Various models have been formulated for the sputter yields 
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(ref. 36 to 39) but again, there Is question as to their accuracy at very low 
Ion energies (ref. 40). Bohdansky (ref. 39) has derived an expression for the 
sputtering yield for the threshold region as: 


where: 


S(E) ~ cc s n (E) 




(26) 


« energy Independent function of the mass ratio between the target ( M 2 ) 
and Ion projectile (M-|) 

E th sputtering threshold energy 


S n nuclear stopping cross section 

Askerov and Sena have proposed an emperlcally formulated expression (ref. 41) 

S(E) oc k (E - E t h) 3 (27) 

Beatle has used an exponential function to approximate Askerov and Sena's data 
to calculate the lifetime of Ion thrusters (ref. 25). 


There Is even considerate disagreement as to what the sputtering thres- 
hold energy should be. Experimental values for Hg -> Mo system (Mo screen 
grid and Hg propellant In an Ion thruster) vary from 16 (ref. 41) to 32 V 
(ref. 42). A theoretical value as calculated by an expression formulated by 
Bohdansky Is as high as 74 V. 

Figure 2 shows the comparison of normalized equations (26), (27), and the 
sputtering yield measurements of Askerov and Sena. For a best fit, a thres- 
hold value of 32 V was used with equation (26) and 16 V as measured by Askerov 
and Sena with equation (27). As suggested by Bohdansky, the S n (E) term was 
assumed to be energy Independent In the Ion energy range of Interest. From 
figure 2, It Is apparent that the best agreement with the measured data Is 
obtained by equation (27) with K = 2 . 61 xl 0 -8 and Eth = 16-0 v - Th1s Implies 
that Bohdansky 's model and Askerov and Sena's data are not In complete agree- 
ment. Optical spectroscopy data from Ion thrusters have Indicated that sput- 
tering threshold energies are considerably lower than the 32 V used In 
equation (26). 


For these reasons, equation 27 represented the best available method for 
determining the sputtering yield In this study, even though theoretically It 
does not seem to have much justification. The sputter yield of double charged 
Ions was assumed to be: 


S++(E) = S + (2E) . 

The accuracy of equations (18) and (27) can be evaluated by comparing the 
calculated values with sputtering rate measurements of long term thruster life 
tests. Table III lists the life tests as summarized by Beatle (ref. 2). Com- 
parison of sputtering rates In table III reveals very good agreement between 
Ion thruster life-test results and the calculated rates using equation (18). 
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Ion Thruster Plasma Parameters 


The comparison of the spectroscopic model to measured Intensity data 
requires the definition of the electron temperature and density as well as the 
doubly to singly charged Ion current density ratio In calculating the sputter- 
ing rate as defined In the previous section. Analytical models have been 
developed to calculate these parameters (refs. 44 and 45), however, these cal- 
culations are rather complex. In this study, the plasma parameters were deter- 
mined from either semlemplrlcal equations or strictly empirical expressions 
obtained from a limited amount of probe data available In the literature. Data 
was characterized as either 8- (refs. 46 to 48) or 30-cm (refs. 1 and 49 to 53) 
mercury Ion thruster data. Possible differences between thrusters of the same 
size as grid and baffle configurations or other geometrical variations were 
Ignored In applying the probe data In the literature to the study discussed 
herein. 

Electron temperature and density . - It has been recognized that In an Ion 
thruster more than one electron density distribution Is possible. Electrons 
may exhibit what Is referred to as primary electron distribution In which the 
electrons have energies that of the plasma potential or It Is possible for the 
electrons to have undergone sufficient thermallzatlon collisions resulting In 
a Maxwellian distribution. Most recent probe data (ref. 1) Indicates that 
Maxwellian electron distributions dominate In both thruster sizes, therefore, 
only the Maxwellian electron distribution will be considered. 

The Maxwellian temperature T^ Is related to the Ionization rate fac- 
tor Q$ defined as: 



where e Is the electronic charge, m e Is the electron mass, oj(E) Is the 
Ionization collision cross section and E Is the electron energy. 

Instead of solving equation 28, the following, much simpler equations, as 
obtained from available probe data, will be used for the respective thruster 
sizes: 


<Ty) - 
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(29a) 


(29b) 


30 9.0x10 

The Maxwellian electron number density can be derived for an Ion thruster as: 

J, 


n.. « 
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(30) 
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Where Jg 1s the ^ on & earn current, Vg Is the Ion Boehm velocity, Ag Is the 
area of the grids and <p s Is the transparency of the grids. 

The following equations, of the form of equation 30, were obtained for the 
Maxwellian electron number density from available data for the 8- and 30-cm Ion 
thrusters : 


(n M ) 


<V 


30 


8.73xlO ? (J B ) 


u 8 

/ n \ 

1 .2xl0 11 (J n ) 

8 30 

(cm 3 ) 

/ -D-\ 


(cm 3 / 


(31a) 


(31b) 


where (Jg)g and (Jg)30 are measured In mA, and amps respectively. It Is 
noted that the ratio of the two electron number densities, as calculated by 
equations (31a) and (31b) for a given beam current Is very close to the ratio 
of the two thruster sizes as predicted by equation 30 (assuming the same grid 
transparency). The 30-cm thruster parameters as expressed by equations (29b) 
and (31b) were obtained from data over a narrow thruster operating range. The 
probe data for varying beam current ranged from 1.0 to 2.0 A at a discharge 
voltage of 32 V. The discharge voltage was varied from 29 to 36 V at a beam 
current of 2.0 A (ref. 52). 

It should also be noted that severe plasma Instabilities were observed by 
reference 47 In measuring the probe data of the 8-cm thruster between discharge 
voltages of 30 to 40 V at discharge currents of greater than 1.0 A. The quoted 
expressions, therefore, for the electron temperature and density may not be 
that accurate over the whole operating range of the 8-cm Ion thruster. 

Doubly to single charged Ion current density ratio . - The doubly to singly 
charged Ion current density ratio, J ++ /J + , can be calculated from the ratio 
of the respective number densities: 


++ 


.++ 


v w A w - v b\ 
F + W 


(32) 


where Q* + Is the doubly charged Ionization rate coefficient, V* + 
are the reduced Ion and Boehm velocity of doubly charged Ions, Aw 

F+ + 


and V 


++ 


B 


effective area for Ion loss at the wall, Ag 
formity factor for the density gradient and 


Is the 

Is the beam area, F ++ Is the unl- 
Q Is the Ion-production volume. 
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Again, for simplicity, the ratio was obtained from available probe data In 
the literature: 


| = 1.67(V d - 28.8) (33a) 

8 


ry U • ur| 

= 1.3x10"^ e u (33b) 


where nu Is the propellant mass utilization In percent. It Is noted that 
the ratio was found not to correlate with the same parameter In both thrusters. 

Triply charged Ions have been observed In 30-cm thruster above discharge 
voltages of 32.5 V (ref. 53). Centerline measurements Indicated ratios of 
only about 10 percent at a discharge voltage of 50 V. This effect, however, 
was neglected. 




EXPERIMENTAL APPARATUS 

Two conventional 8- and 30-cm thrusters were used for the optical spec- 
trometer emission measurements. The 30-cm thruster was equivalent to a ''900" 
series EMT thruster described In reference 54. The 8-cm thruster was slightly 
modified from Its baseline design to be able to operate at higher power levels. 
The description of this thruster In more detail Is found In reference 55. The 
30-cm thruster was operated In the 7.6-m diameter by 21.4-m long vacuum facil- 
ity at NASA's Lewis Research Center. Its no load pressure was about 
l.OxlO -7 torr. The 8-cm thruster was operated In the 1.5 by 6.1 m vacuum 
tank, capable of about 5.0xl0 -7 torr pressure at no load conditions. 

A Jarrel-Ash, Model 82-000 series 0.5 m Ebert scanning emission spectro- 
meter with an electric drive was used with a strip chart recorder to monitor 
the spectral line Intensities emanating from the Ion thruster. The spectro- 
meter was located at the opposite end of the facility In which the thruster 
was operated. The quartz window through which the spectrometer viewed the 
thruster was In a direct path of the beam current Inside the facility. This 
configuration proved to be the most successful In ensuring that the quartz 
window remained clean from sputter deposition and facility oil residue. A 
slit on the spectrometer, usually set about 75 microns, viewed the entire 
diameter of the thruster (about 5 percent of the beam area). The results of 
measurements obtained with the 30-cm thruster In studying the effect of back- 
ground pressure on the sputtering rate of thruster components can be found In 
references 43 and 1. Whereas, various lines of mercury and molybdenum were 
observed, only the Mol 3789 A line Intensities will be discussed herein. 
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RESULTS AND DISCUSSION 


Comparison Between Model and Measured Intensity 

The relative Intensity ratio at two thruster operating points from equa- 
tion 25 becomes: 


I 

I 


1 

2 


[ V 2 R 1 

tn J ] 1 R 2 


(34) 


where R can be thought of as the total number of sputtered atoms (assuming 
constant sputtered Ion velocity) and nj as the fraction of excited sput- 
tered atoms. 


The simplest condition to evaluate In an operating thruster Is when Ion 
energy Is varied without affecting the plasma parameters, namely changing R 
without changing nj. This condition approximates Wehner's experimental 
study of sputtering at low Ion energies with an optical emission spectometer. 

In an Ion thruster, this condition can only be achieved by electrically Insu- 
lating the screen grid from the discharge chamber and then applying a negative 
bias to the screen grid with respect to the discharge chamber. The Ion energy, 
In this case, was calculated by modifying equation (19) as: 


(E) 3 o - !.29 (V D ) - 9. 38 + v bias (35) 

The comparison of the measured relative Intensity of the Mol 3798 A line 
and the calculated sputtering rate ratio with an applied bias to the 30-cm 
thruster screen grid Is seen In figure 3. The direct sputtering rate measure- 
ments of Askerov and Sena and sputter yields as determined by equation (27) 
were used to calculate the two normalized Intensities. It Is apparent that 
excellent agreement Is achieved with the absolute sputtering yields and satis- 
factory agreement Is obtained with equation (27). The discrepancy may partly 
be due to an overestimation of the sputter yield by equation (27) at Ion 
energies of greater than 80 V as seen In figure 2. The comparison, of course, 
Is based on the assumption that changes In the measured Intensity are the 
direct result of changes of the sputtering of the screen grid only. As later 
will be demonstrated, sputtered molybdenum Intensity can change by other 
mechanisms. The relatively good agreement between the calculated and measured 
Intensities validates the use of equation (27) for the sputter yield as well 
as justifying the assumption of constant sputtered Ion velocity with Ion 
energy discussed previously. 

Once the discharge conditions are changed, the physical conditions Inside 
the plasma are changed and the situation becomes much more complicated as 
compared to the simple case just discussed. The relative population levels of 
the excited sputtered atoms nj must be considered. 

The analytical model formulated for nj appears In the first part of the 
paper and listed In table I. Values of nj were obtained for the selected 
transition (In this case the Mol 3798 A line) by the proper choice of electron 
.temperature and density for a particular thruster operating point. The plasma 
parameters were determined by equations (29), (31), and (33) at the given 
thruster operating parameter. 
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It should be pointed out that the plasma parameters from the thruster 
centerline were utilized, whereas, the spectrometer observed the entire diam- 
eter of the thruster. However, comparison of Intensities measured at the cen- 
ter of the thruster with those measured across the entire diameter Indicate 
that this discrepancy should not be of any significance (ref. 1). 

Figure 4 shows the comparison of the measured and calculated relative 
Intensity of the Mol 3798 A line as a function of discharge voltage of the 
30 cm thruster. The Increased Intensity with discharge voltage Is associated 
mainly with the Increasing sputtering rate with Ion energy. The fraction of 
excited sputtered atoms, th, remains relatively constant with discharge 
voltage. The agreement between the model and the measured Intensities Is good 
considering the many assumptions made In the model for calculating n i , and 
the simplification and extrapolation of the plasma parameters by equations 
(29), (31), and (33). The lower than expected Intensity measurements at the 
low Ion energy range may be partly due to the attenuating effects on sputtering 
due to background gases of the vacuum facility. 

Figure 5 shows the relative Mol Intensity as a function of beam current. 

In this case, the Increase In Intensity with Increasing beam current Is asso- 
ciated with both Increased sputtering rate (proportional to the beam current) 
and the Increase In the fraction of excited spultered atoms, nj. The Increase 
In nj with Increasing beam current Is due to an almost linear Increase of 
electron density with beam current. The electron temperature was assumed to 
be constant. Again, part of the discrepancy at the low beam current range may 
be due to the effects of background gases on the sputtering rate as mentioned 
before. The assumption of constant electron temperature may also not be 
strictly valid. 

The baseline 8-cm thruster has only one propellant Inlet as compared lo 
two of the 30-cm thruster. This restriction of the smaller thruster does not 
always allow Independent thruster parameter variation as Is possible with the 
larger thruster. The beam current In the 8-cm thruster Is Increased by 
Increasing the propellant flow rate through the cathode allowing the discharge 
voltage to reach Its desired value. Figure 6 shows the relative Intensity of 
an 8-cm thruster as a function of discharge voltage. The value of the beam 
current Increased from 106.3 mA at at discharge voltage of 30 V to 211 mA at 
25.4 V as Indicated at each data point. The sputtered molybdenum Intensity Is 
expected to decrease with decreasing discharge voltage, but the decrease Is 
moderated by Increasing beam current and associated Increase In n j • As 
mentioned before, the probe data of the 8-cm thruster between 30 and 40 V of 
discharge voltage was severely Impacted by plasma Instabilities, making an 
accurate determination of the plasma parameter difficult. The large discrep- 
ancy at the low discharge voltages was not, however, mainly due to the lack of 
accurate plasma parameters. It was surmised that the large Increase of Inten- 
sity at the low discharge voltage was mainly due to increased Ion Impingement 
on the accelerator grid. It was noted that as the beam current Increased, 
accelerator Impingement current significantly Increased due to poor perveance 
or misalignment of the grids. The evidence leading to this conclusion Is the 
Increasing accelerator current to beam current ratio as the Intensity Is 
observed to Increase. (Mol Intensity has been observed to Increase suddenly 
when Ion beams have been Intentionally defocused.) 

Figure 7 shows the relative Mol Intensity as a function of beam current 
with the discharge voltage constant (this Is possible only at low discharge 


21 



voltages for the 8-cm thruster). The agreement between theory and measured 
Intensity appears to be very good In this case. 


Optical Spectroscopy as a Technique for Determining Thruster Lifetime 

It appears that the model developed to correlate spectroscopic emission 
Intensity to the sputtering rate of a thruster component holds some promise as 
a technique for rapid determination of thruster lifetime. Satisfactory agree- 
ment was attained between measured emission Intensity and calculated Intensi- 
ties of 8- and 30-cm Hg Ion thrusters, despite less than perfect definition of 
plasma parameters and sputtering yields at very low Ion energies. For the 
Hg -> Mo system considered In this study, sputter yields have been measured 
with Ion energies as low as 40 V. For other propellant - target combinations, 
this Is not the case. For example, Xe -* Mo sputtering yields are not availa- 
ble below 100 V. It was concluded before that theoretical models developed 
for sputtering yields at low energies still appear to be Inadequate. Thus It 
Is for such propellant-target combinations that the optical spectroscopy tech- 
nique may offer the most promising possibilities. The requirements for such 
an application are Improved definition of plasma parameters as a function 
thruster operating parameters. The population levels of excited Mol, nj, as 
a function of electron density and temperature, have been already determined 
In this study. The measurement of the excited Mol relative Intensity, there- 
fore, should result In an estimate of a xenon Ion thruster lifetime by calcu- 
lating the sputtering rate R by equation (34). The relative Intensity 
measurements taken below discharge voltages of 100 V can be converted to abso- 
lute sputtering rates by normalizing the measured Intensity to the measured 
sputtering rates at 100 V and above. 


CONCLUSIONS 

An analytical model has been formulated to correlate the measured emis- 
sion Intensity of excited sputtered atoms In an Ion thruster to the sputtering 
of thruster discharge components such as the screen grid. First, however, the 
thruster discharge plasma was characterized In regard to the classical spec- 
troscopic parameters such as equilibrium conditions and optical depth of the 
plasma. It was established that the plasma Is not In Local Thermal Equilib- 
rium (LTE) and that the plasma Is optically thin (no loss of radiation Inside 
the plasma) . 

The model to calculate population densities of various transitions used a 
kinetic steady-state approach balancing the rates of population and depopula- 
tion of ten low-lying excited states of the sputtered Mol atom with those of 
the ground state. This was done to relate the spectral Intensities of the 
various transitions of Mol to the population densities of the ten various 
states of the atom. Once this Is done, the population density can be related 
to the sputtering rate of the target. Two modes of sputtered atom excitation 
and decay (radiative and colllslonal) were considered. 

Since actual data has not been published of Mol excitation rate con- 
stants, semlemplrlcal expressions were used. A computer was utilized to solve 
the 10 by 10 matrix of rate equations In obtaining the excited sputtered atom 
density as a function of electron temperature and density. 
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The calculated relative Intensity of Mol 3798 A was compared to measured 
relative emission Intensities of the same transition from 8- and 30-cm mercury 
Ion thrusters using semlemplrlcal equations for the electron density, tempera- 
ture and sputtering yields of a Hg -» Mo system. Excellent agreement between 
theory and data was found for the special situation In a thruster where the 
Ion energy Is varied without changing the plasma parameters. (This was accom- 
plished by biasing the screen grid negatively with respect to the discharge 
chamber. ) 

Satisfactory agreement between calculated Intensities and relative meas- 
ured Intensities was obtained when Intensities were measured as a function of 
beam current (at constant discharge voltage) and as a function of discharge 
voltage (at constant beam current), for the 30-cm thruster. Agreement between 
the theory and data was not, however, satisfactory for the Intensity measure- 
ments of 8-cm thruster as a function of discharge voltage (with varying beam 
current). However, most of this discrepancy can be related to the emission 
created by Ion Impingement on the accelerator grid due to poor perveance or 
grid misalignment. Part of the discrepancy may be also due to Insufficiently 
accurate definition of electron density and temperature In the 8-cm thruster. 
This was a result of Instabilities experienced during probe measurements over 
a wide range of thruster operating parameters. 

The technique formulated, however, appears to hold promise for rapid 
evaluation of thruster lifetimes, especially In configurations where sputtering 
yields at low Ion energies are not available In the literature (e.g., Inert gas 
Ion thruster systems). 
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APPENDIX A 


Much of the data required for the KSS model Involves properties of the 
screen grid material, molybdenum. The spectroscopic system under considera- 
tion for molybdenum Is shown In figure A(l). Ten levels, representing low- 
lying excited states, are shown, along with the ground state of molybdenum. 
These levels form five sets of levels, the single 7 S ground state level, a 
pair of 5 P levels, two sets of 7 P levels, each containing three levels and 
a pair of 7 D levels. The two sets of 7 P levels are prefixed with "y" and 
"z", representing the lower and upper set of 7 P states, respectively. 

Table A-l presents some Important data for this system. The first column 
shows the term symbol for each of the eleven levels. The second column gives 
the wavelength for the resonance transition, the transition between this level 
and the ground state, In A. Column 3 presents the frequency of that transi- 
tion, In s _1 . In column 4 the statistical weight of the level Is given; a 
dimensionless number. 

Columns 5 and 6 of Table A-l show the energy values for the 10 transi- 
tions comprising this system, In both eV and ergs. Both of these units are 
useful for certain calculations; both are derived from the wavelength data and 
are Included here for ease of reference. 


Equations (6) and (7) of this report require the absorption oscillator 
strength, f^j, as part of their Input data requirements, while equation (8) 
requires the related parameter, Aj^, the Einstein transition probability. 
Several sources tabulate oscillator strengths and transition probabilities 
(refs. 12, 43, 44, 56, and 57), Including one source which attempts to evalu- 
ate the quality of these values (ref. 12). Several works have recently been 
published on the lifetimes of various molybdenum atom states (refs. 58 and 
59). If only resonance transitions are considered significant, the following 
equation converts lifetime to absorption oscillator strength 


f 


U 


1 .499x10 


-7 £l xi 
9 1 T 


(Al) 


where -r Is the lifetime of the state In ns, X Is the wavelength In A, 
gj and g^ are the dimensionless statistical weights of the upper and 
lower states, respectively and f^i is the absorption oscillator strength, 
which Is also dimensionless (ref. bO). If the absorption oscillator strength 
Is known, the Einstein transition probability can be easily found via 


0.6670x10 16 , 

= X 2 9 j ^ 


(A2) 


where x Is In A and Aji Is In s -1 (ref. 60). Table A-2 presents the 
upper state lifetimes, the calculated absorption oscillator strengths and the 
Einstein transition probabilities. For resonance transitions, the Einstein 
transition probabilities are equivalent to the reciprocals of the lifetimes of 
the upper state; this serves as a check of the results of equations (Al) and 
(A2). 
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TABLE AT. - MOLYBDENUM TERM SYMBOLS, WAVELENGTHS, 
FREQUENCIES, AND STATISTICAL WEIGHTS FOR THE 
ELEVEN LEVEL SYSTEM 3 


Term 

symbol 

a 

V, 

Hz 

g 

■ 

E, 

ergs 

7 s 



7 



b 3 








z 7 P 2 

3903.0 

7.6811xlt 

.14 

5 

3.1769 

5 . 0896x1 C 

12 


3864.1 

7.7584 


7 

3.2089 

5.1409 



3798.3 

7.8928 


9 

3.2645 

5.2299 



3466.8 

8.6475 


5 

3.5767 

5.7300 



3456.4 

8.6735 


7 

3.5874 

5.7473 



3208.8 

9.3428 


5 

3.8642 

6.1907 


» 1p 2 

3194.0 

9.3861 


5 

3.8821 

6.2194 


» 7p 3 

3170.3 

9.4563 


7 

3.9112 

6.2659 


7 °s 

3158.2 

9.4925 


7 

3.9261 

6.2899 


A 

3132.6 

9.5701 1 

' 

9 

3.9582 

6.3412 ' 



Reference 11 













TABLE A2. - MOLYBDENUM UPPER STATE LIFETIMES, ABSORPTION OSCILLATOR 
STRENGTHS, AND TRANSITION PROBABILITIES 


Term 

Upper state 
lifetime, 
ns 


Aj 1 10“ 7 
Hz 

f 



(a) 

(b) 

(a) 

(b) 

(c.d) 

(a) 

(b) 

(c,d) 

% 

22.1 


d 4. 525 



d 0. 05823 



% 

22.3 

— 

4.484 


0.2147 

0.08031 


0.003845 


17.1 

16.3 

5.848 

6.135 

9.089 

0.09536 

0.1001 

.1483 


17.0 

15.8 

5.882 

6.329 

6.530 

.1317 

.1417 

.1462 


15.9 

14.7 

6.289 

6.803 

5.267 

.1749 

.1891 

.1465 

a 2 

— 

6.6 


15.15 

20.54 


.1655 

.2244 

> 7 >3 

— 

7.3 



13.70 

12.21 


.2064 

.1840 


— 

5.6 


17.86 

12.86 


.3378 

.2431 

m 

— 

35.6 


2.809 

4.155 


.03097 

.04582. 

\ 



21.2 


4.717 

4.469 . 


.07052 

.06682 


Reference 58. 
bReference 59. 

Reference 61 . 

d Used for generating tables II a to j. 


26 

















REFERENCES 


1. Poeschel, R.L. and Beattie, J.R., "Primary Electric Propulsion Technology 
Study," NASA CR-159688, 1979. 

2. Beattie, J.R., "Endurance Test of a 30-cm Diameter Engineering Model Ion 
Thruster - Investigation of Thln-Fllm Erosion Monitors for Ion Thrusters," 
NASA CR-T 681 32, 1983. 

3. Rock, B.A. and Parsons, M.L., "Rapid Evaluation of Ion Thruster Lifetime 
Using Optical Emission Spectroscopy," NASA CR-174741, 1984. (Not yet 
published.) 

4. Wilson, R., "The Spectroscopy of Nonthermal Plasma," Journal of 
Quantitative Spectroscopy and Radiative Transfer . Vol. 2, 1962, pp. 

477-490. 

5. McWhlrter, R.W.P., In Plasma Diagnostic Techniques . Academic Press, New 
York, 1965, pp. 201-264. 

6. Grlem, H.R., Plasma Spectroscopy . McGraw-Hill, New York, 1964, pp. 196-198. 

7. Alkemade, C. Th. J., Hoolander, Tj., Snelleman, W., and Zeegers, P.J. Th., 
Metal Vapors In Flames . Pergamon Press, Oxford, 1982, pp. 125-129. 

8. Alkemade, C. Th. J., Hoolander, Tj., Snelleman, W., and Zeegers, P.J. Th., 
Metal Vapors In Flames . Pergamon Press, Oxford, 1982, pp. 153. 

9. Alkemade, C. Th. J., Hoolander, Tj., Snelleman, W., and Zeegers, P.J. Th., 
Metal Vapors In Flames . Pergamon Press, Oxford, 1982, pp. 487. 

10. Unpublished data, obtained from M. Mantenleks, NASA Lewis Research Center, 
Cleveland, OH, 44135. 

11. Mavrodlneanu, R. and Bolteux, H., Flame Spectroscopy . Wiley, New York, 

1965, p. 438. 

12. Parsons, M.L., Smith, B.W., and Bently, G.E., Handbook of Flame 
Spectroscopy , Plenum Press, New York, 1975, pp. 295-334. 

13. Lochte-Holtgreven, W., ed., Plasma Diagnostics , North Holland Publishing 
Co., Amsterdam, 1968, pp. 135-213. 

14. Howatson, A.M., Introduction to Gas Discharges . 2nd ed., Pergamon Press, 
Oxford, 1976, pp. 134-135. 

15. Lovett, R.J., "A rate model of Inductively coupled argon plasma analyte 
spectra," Spectrochlmlca Acta , Vol. 37B (11), 1982, pp. 969-985. 

16. Alkemade, C. Th. J., Hoolander, Tj . , Snelleman, W., and Zeegers, P.J. Th., 
Metal Vapors In Flames . Pergamon Press, Oxford, 1982, pp. 107-109. 

17. Wlneforder, J.D., McGee, W.W., Mansfield, J.M., Parsons, M.L., and 
Zacha, K.E., "Intensity of Thermal Radiation of Metal Spectra In Flame 
Emission Spectrometry," Analytlca Chlmlca Acta . Vol. 36, 1966, pp. 25-41. 


27 



18. Peters, R.R., "Double Ion Production In Mercury Thrusters," NASA 
CR-135019, 1976. 

19. M.L. Parsons, private communication. 

20. Bauer, B. and Bartkey, C.D., "Calculation of Inelastic Electron-Molecule 
Collision Cross Sections by Classical Methods," Journal of Chemical 
Physics . Vol . 42, July 1965, pp. 2466-2476. 

21. Sears, F.W. and Zemansky, M.W., University Physics . 4th ed., 

Addlson-Wesley Pub. Co., Reading, MA, 1970, pp. 118-120. 

22. Collett, C.R., "Thruster Endurance Test," NASA CR-135011, 1976. 

23. Collett, C.R. and Bechtel, R.T., "An Endurance Test of a 900 Series 30-cm 
Engineering Model Ion Thruster," AIAA Paper 76-1020, November 1976. 

24. Mantenleks, M. A. and Rawlln, V.K., "Sputtering Phenomena of Discharge 
Chamber Components In a 30-cm diameter Hg Ion Thruster," AIAA Paper 76-988 
November 1976. 

25. Beattie, J.R., "A Model for Predicting the Wearout Lifetime of the 
LeRC/Hughes 30-cm Mercury Ion Thruster," AIAA Paper 79-2079, October 1979. 

26. Oechsner, H., "Sputtering. Review of Some Recent Experimental and 
Theoretical Aspects," Applied Physics . Vol. 8, March 1975, pp. 185-198. 

27. Wehner, G.K. and Anderson, G.S., In Handbook of Thin Film Technology . 
McGraw-Hill, New York, 1970, pp. 1-38. 

28. Greene, J.E. and Barnett, S.A., "Ion-surface Interactions During Vapor 
Phase Crystal Growth by Sputtering, MBE, and Plasma-Enhanced CVD: 
Applications to Semiconductors," Journal of Vacuum Science Technology . 

Vol. 21, no. 2, 1982, pp. 285-302. 

29. Stuart, R.V., Wehner, G.K., and Anderson, G.S., "Energy Distribution of 
Atoms Sputtered from Polycrystalllne Metals," Journal of Applied Physics . 
Vol. 40, no. 2, 1969, p. 803-812. 

30. Wehner, G.K., "Forces on Ion-Bombardment Electrodes In a Low-Pressure 
Plasma," Journal of Applied Physics . Vol. 31, no. 8, 1960, pp. 1392-1397. 

31. Harper, J.M.E., Cuomo, J.J., and Kaufman, F.R., "Technology and 
applications of Broad-beam Ion Sources In Sputtering -2. Applications," 
Journal of Vacuum Science Technology . Vol. 21, no. 3, 1982, pp. 737-756. 

32. Bernhardt, F., Oechsner, H., and Stumpe, E., "Energy Distributions of 
Neutral Atoms and Molecules Sputtered from Polycrystalllne Silver," 

Nuclear Instruments and Methods . Vol. 132, 1976, p. 329-334. 

33. Andersen, H.H. and Bay, H.L., "A Survey of Sputterlng-Yleld Data for 
Plasma-Wall Interaction Calculations," Journal of Nuclear Materials . 

Vol. 93-94 (B), 1980, pp. 625-633. 

34. Behrlsh, R., ed., Sputtering by Particle Bombardment . (Topics In Applied 
Physics, Vol. 47), Sprlnger-Verlag, Berlin; New York, 1981. 

28 



35. Roth, J., Bohdansky, J., and Ottenberger, "Data on Low Energy Light Ion 
Sputtering," IPP 9/26, Max-Planck-Instltut Fur Plasmaphyslc, 1979. 

36. Harrison, Don E., Jr., "Sputtering Models - A Synoptic View," Radiation 
Effects . Vol . 70(1-4), 1983, pp. 1-64. 

37. Yamamura, Y., Matsunaml, N. t and Itoh, N., "Theoretical Studies on an 
Empirical Formula for Sputtering Yield at Normal Incidence, Ra diation 
Effects . Vol. 71(1-2), 1983, pp. 68-86. 

38. Strydom, H.J. and Grles, W.H., "A Comparison of Three Versions of 
Sigmund's Model of Sputtering Using Experimental Results," Radiation 
Effects Letters Section . Vol. 86(4), 1984, pp. 145-151. 

39. Bohdansky, J., "A Universal Relation for the Sputtering Yield of Monatomic 
Solids at Normal Ion Incidence," Nuclear Instruments and Methods In 
Physics Research . Sect. B, Vol. 230(1), 1984, pp. 587-591. 

40. Zalm, P.C., "A Critique of Semlemplrlcal Formulas for the Sputtering Yield 
near Threshold Energy," Radiation Effects Letters Section . Vol. 86(1), 
1983, pp. 29-34. 

41. Askerov, Sh.G. and Sena, L.A., "Cathode Sputtering of Metals by Slow 
Mercury Ions," Sovlet-Phvslcs-Solld State . Vol. 11, no. 6, December 1969, 
pp. 1288-1293. 

42. Stuart, R.V. and Wehner, G.K., "Sputtering Yields at Very Low Bombarding 
Ion Energies," Journal of Applied Physics . Vol. 33, no. 7, 1962, pp. 
2345-2352. 

43. Rawlln, V.K. and Mantenleks, M.A., "Effect of Facility Background Gases on 
Internal Erosion of the 30-cm Hg Ion Thruster," AIAA Paper 78-665, April 
1978. 

44. Longhurst, G., "Prediction of Plasma Properties In Mercury Ion Thrusters," 
NASA CR-159448, 1978. 

45. Brophy, J.R., "Ion Thruster Performance Model," NASA CR-174810, 1984. 

46. Wessel, F.J. and Williamson, W.S., "8-cm Ion Thruster Characterization," 
NASA CR-1 67887, 1983. 

47. Wilber, P., "Advanced Space Propulsion Thruster Research," NASA CR-165584, 
1981 . 

48. Hyman, J., Jr. and Williamson, W.S., "Ion Chamber Sputtering 
Investigation," Hughes Research Laboratory, Malibu, CA, Contract IS-679, 
January 1977. 

49. Poeschel, R.L., "2.5 kW Advanced Technology Ion Thrusters," NASA 
CR-135076, 1976. 

50. Poeschel, R.L., "High-power and 2.5 kW Advanced Technology Ion Thruster," 
NASA CR-1 35163 , 1977. 


29 



51. Poeschel, R.L., "Extended Performance Solar Electric Propulsion Thrust 
System Study," NASA CR-135281, 1977. 

52. Beattie, J.R., private communication. 

53. Vahrenkamp, R.P., "Measurement of Double Charged Ions In the Beam of a 
30-cm Mercury Bombardment Thruster," AIAA Paper 73-1057, November 1973. 

54. Bechtel, T.R. and Rawlln, V.K., "Performance Documentation of the 
Engineering Model 30-cm Diameter Thruster," AIAA Paper 76-1033, November 
1976. 

55. Mantenleks, M.A., "Extended Performance 8-cm Mercury Ion Thruster," AIAA 
Paper 82-1913, November 1982. 

56. Moore, C.E., Atomic Energy Levels . 1st ed., National Bureau of Standards, 
Washington, DC, 1958, Vol . Ill, pp. 1-6. 

57. Corliss, C.H. and Bozman, W.R., Experimental Transition Probabilities for 
Spectral Lines of Seventy Elements , Is ed., National Bureau of Standards, 
Washington, DC, 1962, pp. 204-219. 

58. Duquette, D.W., Salih, S., and Lawler, J.E., "Radiative lifetimes In Mol 
Using a Novel Atomic Beam Source," Physics Letters . Vol. 83A, no. 5, 1981, 
pp. 214-216. 

59. Kwlatkowskl, M., Mlcall, G., Werner, K., and Zimmerman, P., "Radiative 
Lifetimes and Oscillation Strengths In Mol," Physics Letters . Vol. 85A, 
no. 5, 1981, pp. 273-274. 

60. Mavrodlneanu, R. and Bolteux, H., Flame Spectroscopy . Wiley, New York, 
1965, pp. 514-516. 

61. Schnehage, S.E., Danzmann, K., Kuennemeyer, R., and Kock, M., "Oscillator 
Strengths of Neutral and Singly Ionized Molybdenum," Journal of 
Quantitative Spectroscopy and Radiative Transfer . Vol. 29, no. 6, 1983, 
pp. 507-516. 

62. Dulgeroff, C.R., Beattie, J.R., Poeschel, R.L., and Hyman, J., Jr., "Ion 
Thruster System (8-cm) Cyclic Endurance Test," NASA CR-174745, 1984. 


30 



TABLE I. - ABSORPTIVITIES (k max VALUES) 
[( n l)Mo = 2.0xl0 5 atoms/cm 3 .] 


Upper state 

kmax at 500 k, 
0.0431 eV 

kmax at 11600 K, 
1.0 eV 

k max at 139 200 K, 
12.0 eV 

z 7 P 2 Mo I 

1.2xl0 -2 cm-1 

2.4xl0 -3 cm-1 

7.0x10-4 cm-1 

7 

-2 

-3 

_4 

2 P 3 

1.2x10 

2.4x10 

7.0x10 

2?P 4 

1.2X10" 2 

2.5x1 0~ 3 

7.1xl0" 4 

% 

5 . 1 xl 0 — 3 

1 .lxlO -3 

3.1 xlO -4 

5p 3 

3.4xl0" 4 

7.1xlO" 5 

2.0xl0" 5 

S 

4.4X10" 3 

9.1 xlO -4 

2.6xl0 -4 

y7p 2 

2.1 xl 0 -2 

4.5X10" 3 

1 .3xl0" 3 

y7p 3 

1 .8xl0 -2 

3.7xl0' 3 

l.lxlO" 3 

\ 

6.5X10" 3 

1 .3xl0 -3 

3.9xl0' 4 

y7p 4 

2.4X10 -2 

4.9X10” 3 

1.4xl0 -3 

3 P 1 Hg I 

+2 

4.1x10 

8.5xl0 +1 

2.5X10 1 









TABLE II. - RELATIVE POPULATION VALUES, nj 
(a) Electron temperature = 2.0 eV 




log n e 
10.50 

B 

1 

log n e 
11.00 

log n e 
11.25 

log n e 
11.50 

1 .00e-0 

1 .00e- 

-0 

1 .00e-0 

1.00e-0 

1 .00e-0 

1 .45e-6 

2 . 58e- 

-6 

4.59e-6 

8.16e-6 

1 .45e-5 

2.00e-6 

3.55 


6.32e-6 

1.12e-5 

2.00e-5 

2.74e-6 

4.86 


8.65e-6 

1.54e-5 

2.73e-5 

3.63e-6 

6.46 


1 .1 5e-5 

2.04e-5 

3.63e-5 

4.26e-6 

7.58 


1 .35e-5 

2.40e-5 

4.26e-5 

9.23e-7 

1 .64 


2.92e- 6 

5.19e-6 

9.22e-6 

1 .23e-6 

2.20 


3.92e-6 

6.97e-6 

1 .24e-5 

1 .49e-6 

2.65 


4.71e-6 

8.37e-6 

1 .49e-5 

9 Ale-1 

1 .68 


3.00e-6 

5.32e-6 

9.47e-6 

1 .21 e-6 

2.16 1 

r 

3.83e-6 

6.82e-6 

1 .21e-5 



(b) Electron temperature = 3.0 eV 


1.00e-0 1.00e-0 1 .00e-0 1.00e-0 1.00e-0 9.99e-l 

2.81e-6 5.00e-6 8.89e-6 1.58e-5 2.81e-5 4.99e-5 

3.88 6.90e-6 1.23e-5 2.18 3.88 6.89e-5 

4.96 8.82e-6 1.57e-5 2.79 4.96 8.81e-5 

6.62 1.18e-5 2.09e-5 3.72 6.62 1.18e-4 

7.85 1.40e-5 2.48e-5 4.41 7.84 1.39e-4 

1.88 3.34e-6 5.95e-6 1.06 1.88 3.34e-5 


3.34e-6 5.95e-6 1.06 1.88 3.34e-5 
4.51e-6 8.02e-6 1.43 2.54 4.51e-5 
5.46e-6 9.71 e-6 1.73 3.07 5.46e-5 
3.42e-6 6.08e-6 1.08 1.92 3.42e-5 
4.43e-6 7.87e-6 l'40 2.49 T 4.42e-5 


9.99e-l 
4. 58e-5 
6.31e-5 
8.64e-5 
1 .1 5e- 4 
1 .35e- 4 
2.92e 5 
3.92e-5 
4.70e- 5 
2.99e-5 
3.83e-5 


9.99e-l 
8.88e-5 
1 .22e-4 
1 .56e-4 
2.09e-4 
2.48e-4 
5.94e-5 
8. Ole-5 
9.70e-5 
6.08e~5 
7.86e-5 






















TABLE II. - Continued. 


(c) Electron temperature = 4.0 eV 



log n e log n e log n e log n e log n e 

10.50 10.75 11.00 11.25 11.50 


EHI 


og n e 
2.00 


1.00e-0 1.00e-0 1.00e-0 1.00e-0 l.OOe-O 9.99e-l 9.98e-l 

3.97e-6 7 .06e-6 1.26e-5 2.23e-5 3.97e-5 7.05e-5 1 . 25e-4 


5.48e-6 9.75e-6 1.73e-6 3.08 
6.77e-6 1.20e-5 2.14e-5 3.81 

9.07e-6 1 .61 e-5 2.87e-5 5.10 

1.08e-5 1 .92e-5 3.41e-5 6.07 

2.72e-6 4.84e-6 8.61e-6 1.53 

3.68e-6 6.55e-6 1.1 6e— 5 2.07 

4.48e-6 7.96e-6 1.42e-5 2.52 


5.48e-5 9.74e-5 1 . 73e-4 


6.77e-5 1 .20e-4 2.14e-4 


9.06e-5 1 .61e-4 2.86e-4 


1 .08e-4 1 .92e-4 3.41e-4 


2.72e-5 4.84e-5 8.60e-5 


3.68e-5 6. 54e-5 1 .16e-4 


4.48e-5 7.95e-5 1 . 41 e-4 


2.78e-6 4.95e-6 8.80e-6 1.56 2.78e-5 4.94e-5 8.78e-5 


3.62e-6 6.44e-6 1.14e-5 2.03 


3.62e-5 6.43e-5 1 .14e-4 


(d) Electron temperature = 5.0 eV 




























TABLE II. - Continued. 


(e) Electron temperature = 6.0 


Term 

log n e 
10.50 

m 

log n e 
11.00 


log n e 
11.50 

m 

log n e 
12.00 

% 

1 .00e-0 

1 .00e-0 

1 .00e- 

-0 

l.OOe- 

0 

9.99e-l 

9.99e-l 

9.98e- 

1 

% 

5.72e-6 

1 .02e-5 

1 .81e- 

-5 

3.21 e- 

5 

5.71e-5 

1 .02e-4 

1 .80e- 

4 

% 

7.90e-6 

1 .40e-5 

2.50 


4.44 


7.90e-5 

1 .40e-4 

2.49 



9 . 44e-6 

1 .68e-5 

2.98 


5.30 


9.43e-5 

1 .68e-4 

2.98 



1 .27e-5 

2.25e-5 

4.01 


7.12 


1 .27e-4 

2.25e-4 

3.99 



1 .52e-5 

2.70e-5 

4.79 


8.52 


1 . 51 e-4 

2.69e-4 

4.78 


v 7 p 2 

4.02e-6 

7.16e-6 

1 .27 


2.26 


4.02e-5 

7.15e-5 

1.27 


A 

5.46e-6 

9.70e-6 

1.72 


3.07 


5.45e-5 

9.69e-5 

1.72 



6.66e-6 

1 .18e-5 

2.10 


3.74 


6.65e-5 

1.18e-4 

2.10 



4.11 e-6 

7.30e-6 

1.30 


2.31 


4.10e-5 

7.29e-5 

1.30 



5.37e-6 

9 . 54e-6 

1.70 ' 


3.02 ' 


5.36e-5 

9.53e-5 

1.69 ' 



(f) Electron temperature =7.0 


1.00e-0 

1 .00e-0 

1 .00e- 

-0 

1 .00e- 

0 

9.99e-l 

9.98e-l 

9.97x10-1 

6.39e-6 

1 .14e-5 

2.02e- 

-5 

3.59e- 

5 

6.38e-5 

1 .13e-4 

2.01x10-4 

8.83e-6 

1 .57e-5 

2.79 


4.96 


8.82e-5 

1 .57e-4 

2.78e- 

4 

1 .04e-5 

1 .86e-5 

3.30 


5.87 


1 .04e-4 

1 .85e-4 

3.29 


1 .40e-5 

2.50e-5 

4.44 


7.88 


1.40e-4 

2.49e-4 

4.42 


1.68e-5 

2.99e-5 

5.31 


9.44 


1 .68e-4 

2.98e-4 

5.30 


4 . 53e-6 

8.05e-6 

1 .43 


2.54 


4.52e-5 

8.04e-5 

1 .43 


6.14e-6 

1.09e-5 

1.94 


3.45 


6.14e-5 

1.09e-4 

1 .94 


7.50e-6 

1 .33e-5 

2.37 


4.22 


7.50e-5 

1 .33e-4 

2.37 


4.62e-6 

8.21 e-6 

1 .46 


2.60 


4.61e-5 

8.20e-5 

1 .46 


6.04e-6 

1 .08e-5 

1.91 ’ 

' 

3.40 1 


6.04e-5 

1 .07e-4 

1.91 ' 

' 






















TABLE II. - Continued. 


(g) Electron temperature = 



log n e log n e log n e log n e 

10.50 10.75 11.00 11.25 




g n e log n e 
.75 12.00 


1.00e-0 l.OOe-O 1.00e-0 1.00e-0 9.99e-l 9.98e-l 9.97e-l 

6.96e-6 1 . 24e-5 2.20e-5 3.91e-5 6.95e-5 1.24e-4 2.19e-4 

9.62e-6 1 .71 e-5 3.04 5.41e-5 9.61e-5 1 . 71 e-4 3.03 

1 .1 3e-5 2 .01e-5 3.57 6.35e-5 1.13e-4 2.01e-4 3.56 

1.52e-5 2.70e-5 4.80 8.54e-5 1.52e-4 2.70e-4 4.79 

1 . 82e-5 3.24e-5 5.76 1 .02e-4 1 .82e-4 3.23e-4 5.74 

4.96e-6 8.82e-6 1.57 2.79e-5 4.96e-5 8.81e-5 1.56 

6.74e-6 1 . 20e-5 2.13 3.79e-5 6.73e-5 1.20e-4 2.12 

8.23e-6 1.46e-5 2.60 4.63e-5 8.23e-5 1.46e-4 2.60 

5 .06e-6 8.99e-6 1.60 2.84e-5 5.05e-5 8.98e-5 1.59 

6.63e-6 1 ,18e-5 2.10 3.72e-5 6.62e-5 1.18e-4 2.09 ,F 


(h) Electron temperature = 9.0 


1.00e-0 1.00e-0 l.OOe-O 9.99e-l 9.99e-l 9.98e-l 9.97e-l 

7.46e-6 1.33e-5 2.36e-5 4.19e-5 7.45e-5 1.32e-4 2.35e-4 

1.03e-5 1.83e-5 3.26 5.80e-5 1.03e-4 1.83e-4 3.25 


1 .20e-5 2 . 1 4e-5 3.81 


6.77e-5 1 .20e-4 2.14e-4 3.79 


1 ,62e-5 2.88e-5 5.12 9.10e-5 1.62e-4 2.87e-4 5.10 

1.94e-5 3.46e-5 6.14 1.09e-4 1 .94e-4 3.45e-4 6.12 

5.34e-6 9.49e-6 1.69 3.00e-5 5.33e-5 9.48e-5 1.68 

7 . 25e-6 1.29e-5 2.29 4.08e-5 7.24e-5 1.29e-4 2.28 

8.87e-6 1 . 58e-5 2.80 4.98e-5 8.86e-5 1.57e-4 2.80 

5.44e-6 9.68e-6 1.72 3.06e-5 5.44e-5 9.66e-5 1.72 


7 .14e-6 1.27e-5 2.26 t 4.01e-5 7.13e-5 1.27e-4 2.25 


























TABLE II. - Concluded. 


(1) Electron temperature = 10.0 


1 .71e-5 



log n e log n e log n e log n e 

10.75 11.00 11.25 11.50 


1.00e-0 1.00e-0 1.00e-0 9.99e-l 9.99e-l 
7.90e-6 1.40e-5 2.50e-5 4.44e-5 7.89e-5 


1.09e-5 1.94 3.45 


1 .27e-5 


6.14e-5 1.09e-4 
7.14e-5 1.27e-4 

9.60e-5 1.71e-4 


2.05e-5 3.65 6.48 1.15e-4 2.05e-4 


5.67e-6 1.01 1.79 
7.71 e-6 1.37 2.44 
9.43e-6 1.68 2.98 
5.78e-6 1.03 1.83 


7 . 59e-6 


3.19e-5 5.67e-5 
4.33e-5 7.70e-5 
5.30e-5 9.42e-5 
3.25e-5 5.78e-5 
4.26e-5 7 . 58e-5 



(j) Electron temperature = 11.0 


























TABLE III. - CALCULATED AND MEASURED SCREEN GRID SPUTTERING RATES 3 


Date 

Location 

Size, 

cm 

Duration, 

hr 

Beam 

current, 

A 

Discharge 

voltage, 

V 

Measured 

screen 

erosion, 

A/hr 

Uncertainty, 

A/hr 

Calculated 
screen grid 
erosion. 
A/hr 


Hugl 

les 

3( 


10 000 

1.4 

37.0 






1976 





4 165 

2.0 

36.1 

312 

-- 

295 

mmm 





590 

2.0 

36 



295 

KHMM 


' 



937 

2.0 

32 



80 

1981 

XEOS 



4 263 

1 .96 

32 

64 

- 

80 


Hughes 13 

8 

9 489 

.078 

36.6 


-• 

79 


a As compiled by reference 2. 
^Reference 62. 





















Figure 1. - A Grotrian diagram for a hypothetical 4 level 
system, showing the resonance transitions. Rate con 
stants for the transitions are shown as well. 




Figure 2. - Comparison of Hg— Mo calculated and measured 
sputtering yields as a function of ion energy. 



Figure 3. - Relative Mol intensity as a function of applied 
bias to the screen grid (Jg ■ 2.0 A, AVj ■ 50.0 V). 



RELATIVE Mol (3798 A) INTENSITY 



Figure 4 - Relative Mol intensity as a function 
of discharge voltage (30 cm thruster, J B • 3.0 A). 


RELATIVE Mol (3798 A) INTENSITY 



Figure 5. - Relative Mol intensity as a function of 
beam current, Jg (30-cm thruster, AV ■ 50 V). 


RELATIVE Mol (3798 A) INTENSITY 



DISCHARGE VOLTAGE, V 

Figure 6. - Relative Mol intensity as a function of discharge 
voltage (8-cm thruster, discharge current 3. 0 A. 
designates accelerator to beam current ratio in 
mA. 


RELATIVE Mol (3798 A) INTENSITY 



Figure 7. - Relative Mol intensity as function 
of beam current (8-cm thruster, AVj-26. 0 V). 



GROTRIAN DIAGRAM OF MOLYBDENUM 
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Figure Al. - Grotrian diagram of molybdenum for the system 
under consideration, shewing the ten resonance transitions 
for the ten excited states chosen for use in the model. 
Adapted from reference 11. 
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